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ABSTRACT 

We use a large N renormalization group (RG) method to study a model of 
interacting boson system with a quenched random potential. In the absence of im- 
purities, the pure boson system has a critical point that describes the superfluid- 
Mott-insulator (SF-MI) transition. The SF-MI transition of d dimensional bosons 
belongs to the [d -f 1) dimensional XY model universality class. In this paper, 
we study the dirty-boson critical points in the neighborhood of this pure SF-MI 
critical point. In general, the on-site random potential in the original lattice 
model gives two types of randomness in the effective field theoretic action. One 
is the randomness in the effective on-site repulsion tx;(x) and the other is the 
randomness of the chemical potential w(x). It turns out that d = 2 is the critical 
dimension for both types of disorder but the roles of these two types of disoder are 
reversed as (i = 2 is crossed. Applying e = d — 2 expansion, we found coupled RG 
equations for both kinds of randomness which reveal several non-trivial critical 
points. All the weak random fixed points we found have three or more relevant 
directions. We conclude that the direct SF-MI transition is unlikely to occur near 
two dimensions. 

PACS numbers: 67.40.Db, 05.30.jp, 74.80.Bj 
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I. INTRODUCTION 

The problem of repulsively interacting bosons in a random potential has been the 
subject of intense research recently [1-10]. This so-called dirty boson problem contains the 
essential difficulty of understanding interplay between interaction effect and randomness. 
One of the reasons why this problem is very challenging is that there is no sensible non- 
interacting limit for disordered bosons. That is, the zero interaction limit is pathological in 
the sense that the bosons will condense into the lowest localized state around a small region. 
For the metal-insulator transition of electrons, which has been understood better, disorder 
alone or interaction alone can induce localization of electrons and drive the electrons to 
the Anderson insulator (AI) or the Mott insultor (MI) state [11]. The interplay between 
these two effects has been studied during the last decade but there are still open problems 
[11]. 

The dirty boson problem has direct experimental realization on ^ife in Vycor glass 
and in other porous media [12,13]. This disordered interacting boson problem may be 
used to understand the superconductor-insulator transition in the disordered thin films 
[14] and short coherence length superconductors [15]. Recently Wen and Wu [16] showed 
that the superfluid-Mott insulator transition of bosons with the Chern-Simons gauge field 
can describe the transition between quantum Hall (QH) states in the absence of disorder. 
Therefore, the dirty boson problem with the Chern-Simons gauge field is intimately to the 
quantum Hall transitions in the presence of random impurities. The QH-MI transition for 
the pure system is also studied by Chen et al. [17] in which they studied fermions with 
the Chern-Simons gauge field. 

Following Fisher et al. [1,2], we can write the Hamiltonian of the interacting lattice 
bosons in a random on-site potential as 

H = Hq -\- Hi , 
Hq = - ^ (-Jo A* Sjii) hi ^- '^V ^ hi{hi-l) , 

where hi = and bj is the boson creation operator at the site i. /j, is the average 
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chemical potential that fixes the boson density and 6fii is the random on-site potential 
with zero average. is the hopping matrix element and Jq = J2j Jij ■ order to study 
the critical phenamena of the system, it is convenient to find the efi^ective field theory. We 
will summarize the approach of Fisher et al. [1,2]. First the off-site hoping term in Hi is 
decoupled by introducing the Hubbard-Stratanovich field then the resulting action can 
be expanded in terms of '^i. Since '^i is linearly related to {hi) for small the field '^i 
can be identified as a superfluid (SF) order parameter. It was shown [2] that the effective 
field theoretic action is given by 

' (2) 

i i 

where Wi, Ui are random functions of i with zero average. 

It is established [2] that SF-MI transition of the pure system has two uninversality 
classes, the commensurate case and the incommensurate case. In the commensurate case, 
the superfluid density commensurates with a periodic potential. In this case, SF-MI tran- 
sition is described by a tricritical point which belongs to the universality class of {d + 1) 
dimensional XY model with the dynamical exponent z = 1. In the incommensurate case, 
the SF-MI transition happens on a line in the w — u plane. It is argued [2] that the generic 
SF-MI transition should be the latter case ( with 2; = 2 at the transition ) rather than the 
former case. 

One natural question is that how disorder affects these two different SF-MI transitions. 
The destruction of superfluid in the presence of disorder brought the concept of the Bose 
glass (BG) phase [4,18] in which bosons are localized by disorder. In seminal papers. 
Fisher and coworkers [1,2] suggested a scaling theory of SF-BG transition. They argued 
that superfluid-insulator tansition should occur through the Bose glass phase and the 
generic transition should be described by the action (2) with u + Ui ^ which does not 
have space-time isotropy. In the scaling theory, they postulated that the compressibility 
is totally due to the phonon mode and one of the main result of this assumption is that 
the dynamical exponent z = d. The simulation of the quantum rotor model [3] and some 
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renormalization group calculation [7] partially supported this picture although a recent 
quantum Monte-Carlo calculation [9] contradicts to these results. The earlier work of Ma, 
Halperin and Lee (MHL) [4] was reexamined and the importance of the term that is linear 
in u3 is emphasized. However, it was also argued that MHL theory may apply to the 
possible direct SF-MI transition at the commensurate case ( particle-hole symmetric case 

)• 

We can see that the dirty boson effective action has the strict particle-hole symmetry 
if u -|- = 0. It was also argued that the general commensurate case corresponds to the 
weaker particle- hole symmetric case [2], i.e., u = but Ui ^ 0. Are the transitions of 
particle-hole symmetric and asymmetric cases in the same universality class ? Originally 
Fisher and coworkers [1,2] prefered that even arbitrarily weak disorder will induce the 
Bose glass phase for both of the incommensurate and commensurate cases (Fig.l (a)). 
However, numerical calculations in [8] have not revealed the intervening Bose glass phase in 
the commensurate case, although the superfliud-insulator transition indeed occur through 
the Bose glass phase in the incommensurate case (Fig.l (b)). Singh and Rokhsar [5] 
performed a real space renormalization group analysis for the commensurate case and 
they found the direct transition from SF to MI when the disorder is sufficiently weak 
; the Bose glass is found beyond a threshold ( Fig.l (a) or Fig.l (b) depending on the 
impurity strength). Zhang and Ma [6] considered hard core bosons with disorder. In this 
real space renormalization group analysis of a quantum spin- 1 XY model with transverse 
random field ( the hard core boson model is mapped to this model), they concluded that 
commensurate and incommensurate cases are in the same universality class and the SF- 
insulator transition occurs alwalys from the Bose glass phase (Fig.l (a)). The conventionl 
renormalization group calculation by Weichman and Kim [7] in which they used double 
dimensional expansion around d = 4 partially supported the original picture [1,2] for the 
general dirty boson problem although some technical problems exist. 

In this paper, we are going to study a large N generalization of the original action. 
By doing expansion, we can treat the interaction non-pertubatively in the coupling 
constant. Both types of disorder are assumed to be weak and we do the pertubation in 
the strength of two types of disorder. The large N generalized action of the original dirty 
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boson model in the Eucledian space is given by 



+ J d'^x dr [ (« + «(x)) 4dr4>i +w{^) 4>l4>i + ^ {(pUif 



(3) 



where i = 1,---,N and N = 1 corresponds to the original model. tt(x) and w{'x.) 
are gaussian random functions of x with zero mean and their variances are given by 
(u(x)w(y)) = Uq 5'^(x — y) and (w(x)w(y)) = Wq 5'^(x — y) respectively, where (■ ■ ■) 
means the random average. 

When u+u{x) = iD+w{'x.) = 0, (3) describes the muticritical point of SF-MI transition 
of the pure boson system. Near d = 2, both of u and terms are strongly relevant. In 
this paper, we will take u = w = and study the RG flow of C/q, Wq and vq. The physical 
meaning of setting w = is the following. We tune the chemical potential ^ to make the 
average boson density alwalys commensurate with the lattice. Thus, we will call (3) with 
u = the commensurate dirty boson theory. 

Our RG calculations are done at the critical point and with the renormalized mass term 
w = in the course of the renormalization. Since the effect of 0^ term is calculated exactly 
at each order of 1/N, the coupling constant go is not renormalized. Following Rcf. [19] 
and using dimensional regularization, we move an infrared scale in the renormalized 
theory with fixed bare parameters to obtain the RG flow of the renormalized parameters. 
Introducing dimensionless measures of disorder, U = and W = we 

performed e = d — 2 expansion. The resulting renormalization group equation up to (;^)° 
order is found to be 

dW ~ ~ ~ ~ , 
^ eW + aUW + bW^ , 

dl 

^ = S + bWU - atj^ , (4) 
dl 

= -bU + cW, 

dl 

where a = tt-, 6 = — , c = — and / is the logarithmic measure of the RG flow. This is the 
central result of this paper. 

Looking at e > case, we can immediately see that there is only the pure fixed point 
which is given by fj* = 0, W* = 0. At this trivial flxed point, U is irrelevant and W is 



relevant so that the RG flow goes to the strong randomness regime where our RG breaks 
down. However, for e < 0, there are three fixed points. AU of them have at least one 
relevant direction in the U — W plane. Thus, including u and w, these fixed points in 
the original theory (3) has at least three relevant directions. Therefore, in both cases, the 
direct SF-MI transition in the W — U plane is unlikely to occur due to the absence of 
the weak random fixed point with two or less relevant directions. The superfluid insulator 
transition is alwalys governed by a strong random fixed point which cannot be reached by 
weak randomness expansion. More details of the RG flow will be discussed later. 

The organization of this paper is as follows. In section II, we consider a restricted 
model in which ■u(x) = and show the basic formalism we used. Here, we also examine 
some possible effects of long range interactions. The RG calculation for the generic com- 
mensurate dirty boson problem up to order is presented and the RG equation is 
calculated in section III. we also discuss about the results. In section IV, correction 
due to (f)^ interaction is considered. In section V, we summarize and conclude this paper. 

II. RENORMALIZATION GROUP ANALYSIS FOR 
THE STRONG PARTICLE-HOLE SYMMETRIC MODEL 

In this section, we consider a rather restricted model in which we set 'u(x) = 0. This 
corresponds to the strong particle- hole symmetric case in the sense that more general model 
requires the particle-hole symmetry only in the average sense, i.e., (tt(x)) = and allows 
the local breaking of the symmetry at each site. We would like to consider more general 
long range interactions / d'^x d'^y (/)^(x) V{x—y) 4>^{y) with F(x) oc (7o/|x|''~'^ and V{q) = 
go/q^ ( < A < 1/2 ). The A = case corresponds to the usual short range interaction. 
However, it should be mentioned that these are not the true long range interactions of 
bosons because the true one should be F(x — x')p(x)p(x') and p(x) = i{(f)^do(f) — 9o0V)- 
But the long range interaction between densities will induce the long range interaction in 
the 0^ term. 

In order to perform the RG calculation in the critical theory, we need to introduce 
an arbitrary infrared mass parameter fj, [19]. We will use the dimensional regularization 
method to calculate relevant divergent diagrams. The renormalization of the theory is 
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given by the renormalization of two point T^'^^ and four point F^^) vertices. We will take 
the inverse of the full propagator as a two point vertex and the two-boson scattering 
amplitude as a four point vertex function. The relation between the bare theory and the 
renormalized theory is given by 

rJ,arU^>^;A) = Z-''/\^/^^) rW(g,a;;/.) , (5) 

where T^^^ and T^^^ represent the bare and the renormalized vertices respectively. We 
found that appropriate renormalization condition for F*^^^ can be chosen as 

d 

-n-' = 1 , 

= . 

q=H, lima_>o u>=avn 



d 



r(2) 
r(2) 



Renormalization condition for the scattering amplitude will be discussed later. Also, fol- 
lowing standard procedure, we require the independence of the bare theory with respect 

to jJL, 



d 

^ djji 



4::i,=o, (7) 



A 

where A is the mass parameter of the bare theory. 

Let us start with the evaluation of the self energy to the (1/A^)° order. The {l/Nf 
order self energy diagram is given by Fig. 2 (a). The polarization bubble ni(g, a; = 0) in 
Fig. 2 (b) is calculated as 

_ /• rf'^A; rfz/ 1 1 

3^ ' 



^0 



_ 1 n^)TH^) 

Assuming 3 -|- A > cZ, the diagram in Fig. 2 (c) can be approximated as 

^0 (T^WT-^-^^^r^Y " W^o^g^(^+"-^) . (9) 
\1 + Ui{q,u} = 0){go/q>^)J cfg^ 
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Using this result, we can evaluate the self energy as 



C2 



C3 



1 r(d/2 - 2 - A) r(3 + A - d/2) V{d/2 - 1) 



(47r)^/2 r(rf-3- A) r(2 + A) 

2 + A - d/2 r(d/2 - 2 - A) r(3 + A - d/2) T{d/2 - 2) 



(10) 



(47r)'^/2 r(d - 3 - A) r(l + A) 

where voq ui is assumed. The bare two point vertex T^^^.^ up to (;^)° order is 

r£l=^' + ^oV-Si(?,a;) 

_ 2/-, Wpyj C3 2+2A-d^ 

+ .oV(l-^V+^^-'^). 



(11) 



Note that C2 and C3 diverge at d = 2 + 2A. Therefore, let us try e = d — 2 — 2A expansion 
in order to handle these divergences. Also, for convenience, let us introduce dimensionless 
measures of the disorder Wq = ^2+2\-d ^^le bare theory and W — in 
the renormalized theory. Adding appropriate counter terms to cancel the j divergences 
from C2 and C3 and using the renormalization conditions (5) (6), we get the following 
equations 

cf 11 
vl ^ Z-\^ ( 1 _ ^ In 

where £2 = ec2 and £3 = ec^. From jj, independence of the bare parameters, we get 



^^(InZ) ^ -W% 



C2 _ C3 
^2 „2 



(13) 



Wv . 



Now we are going to renormalize the four point function T^'^\ We will take the 
scattering amplitude of two bosons at qi = = q, q2 = qa = 0,u!i = UJ2 = u!3 = ca^ = u! 



as our r*^'^). First of all, let us calculate F^^'^.g in the bare theory, r^^^-e (;^)° 
order is the sum of T^^L,i^T^hL,2^ ^hL,3^ ^hL^ which are shown in Fig.3 (a), (b), (c), (d) 
respectively. Let us introduce a renormalized polarization buble Hi, a new polarization 
bubble 112 and a vertex Vi which are given by the diagrams in Fig.3 (e), (f) and (g). F^^g^j.^ 
can be calculated symbolically as 



p(4) ^^(4) , ^(4) ^(4) p(4) 

bare bare,l bare, 2 bare, 3 bare, 4 



1 / 1 



, Wo + 2Vi ^ Wo . 

,i + ^o(ni + n2); Vi+^onJ 

The bubble Hi which is renormalized by the self energy correction is given by 

, . d'^k diy Z Z 
ni(g,a; = 0) 



Z'^q^-' 



Evaluation of 112 (g, a; = 0) is a little bit long task and the result is 
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where 

1 F(l/2-A) 



X 



X 



(47r)d+i/2 r(rf-3-A) 

dx2 4/'-'-^ (1 - x^f-^-"^^ 

Jo Jo 

/ dyi dy2y2^~^~^ {x2 + y2{'^-X2))~^''^ 

Jo Jo 



(14) 



(15) 



glcl J {2txY i^TxY 27r vl{k - qf + v^k^ + 1^ 
1 1 1 

^ vl{k-pY + v-^ vl{k - p - qf + p^id-3-X) '^^^^ 

- ^0 ^ „d-3„2+2A-d 
2 9 y y : 



X [xi{l - xi){l - X2)y2 - xlyjx2 - 2yiy2XiX2{l - X2) + 2/1(1 - 2/1)^2(1 - X2)'^] 

(17) 

We found that A ^ and A = cases should be treated separately. First, for A = case, it 
is found that the constant e does not diverge as e ^ so that it can be dropped in the final 
RG equations. Now, let us look at the case of A 7^ 0. The strategy is that we investigate 
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the most divergent contributions in various limits and add up all of the contributions. 
Since the most divergent contribution comes from 0:2 — > or ?/2 ^ limit in Eq.(17) and 
we want just this contribution, we can set 0:2 = or 1/2 = inside the square bracket in 
Eq.(17). We found that both ways give the same answer. Here we will take X2 = limit 
inside the square bracket and multiply by 2. The equation (17) becomes 

2 r(l/2-A) 



e = 



X 



(47r)d+i/2 r{d-3-X) 

)^-i/2 f'~^' dx2 xi/^-^-\l - x^f-^-"/^ 
Jo 



dxi Xi ^^^(1 — xij 



X 



dyi 
2 



dy2 V2 



d/2-3 



r(i/2-A) 



(18) 



X 



X 



(47r)rf+V2 r(d - 3 - A) {d/2 - 1 - A) {d/2 - 2) 

t dx^ xt-^/^l-xi)'^/2-3/2 
Jo 

f dy, (1 - yrY'^-' , 
^0 



where only the leading divergence is taken in the second equation. The evaluation of the 
remaining integrals is straightforward and the result is 

2 r(l/2 - A) V{d/2 - 1/2) r(A + 1/2) V{d/2 - 1) 



- (47r)d+i/2 (^/2 - 1 - A) {d/2 - 2) r(d - 3 - A) V{d/2) V{\ + d/2) 
_2 2 r(l/2-A) r(l/2 + A + e/2) r(A + l/2) r(-l + A + e/2) 

- 7 (47r)d+i/2 r(-l + A + e) r(l + A + e/2) r(l + 2A + e/2) 

Vi(gi = 93 = 9, 92 = 0,a;i = uj^ = U},UJ2 = 0) can be evaluated similarly as 

WqvI f d'^p 1 1 1 



(19) 



V^ = 



'I [ 1 

I J {2^Y vl{q-p 



^0^0 2+2A-d 

,2 y 



dx 



9o 
1 



di 



^1 



pY + VqP'^ + p2(d-3-A) 

2 ^ do 



00 



(20) 



r(d/2 - 1 - A/2) r(2 - d/2 + A) r(d/2 - 1) 



d2 = (iiA 



(47r)'^/2 r(d-2- A) r(l + A) 

1 r(e/2) r(l-e/2) r(A + e/2) 
(47r)^/2 r(A + e)r(l + A) 

2(A-1) d-3-A" 



(21) 



d-4 



d-2- A 
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where voq » a; is assumed. Note that d2 = for the short range interaction (A = 0). 

Putting together all of these results, we can get the following ^Y^^^iqi = = Q2 
qs = 0, cji = CJ2 = 1^3 = (^4 = , 



.(4) 
bare 



1 



^o(ni + n2) 



Wo + 2Vi 



1 



Wn 



WqV^ „2(3-d) 7-4 



Wo e 2+2X-d 

2 3 ^72 y 



+ 2 



Wpyj 
9l 



^2+2A-d ( ^ + 



1 9 



do 



(22) 



The renormalization condition for F^^^ is taken as 



Wv 

^^'^Hqi ^ Q4:^ IJ', ?2 = ?3 = 0, lim ioi ^ U2 = UJ3 ^ C04 ^ av/j,) = ^ ^ . (23) 

We introduce again = ^:^^//^+^'^~'^ and Wq = A^+^^~'^. If appropriate counter 
terms were added to cancel the 1/e divergences, the following equation can be obtained 



Wo ^ WZ' 



d-2-2A 



1 + 2W 4 In ^ - 2W ^ In - 



w 



a; 



d-2-2A 



W e ~ A 
1 + ^ (C3 + C2 + 2 2di) In - 



(24) 



where e = ee and di = ed\ . Using the fact that the bare parameters are fixed, as we change 
fi, we get the following equation for W 



C2 + C3 + 2 2di 

Cl 



(25) 



Let b be the reduction factor for the momentum scale from /i to then the RG equations 
for W and v up to (;^)° order is 



dW 
"df 

dv 

m 



= -P{W) = eW 



C2 + C3 + 2 2di 

Cl 



1 



(26) 
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where I = \nb is the logarithmic measure of the RG flow. £2 and £3 are given by 

2 r(2-e/2) r(A + e/2) 



C2 



C3 = 



(47r)rf/2 r(-l + A + e) r(2 + A) ' 
2-e r(2-e/2) r(-l + A + e/2) 
(47r)'^/2 r(-l + A + e) r(l + A) 



Now, the RG equations for < A < 1/2 up to (;^)° order is 
dW 



= -/3iW) = eW- o ,,,, , ,, W' , 



dl ' c? (47r)^/2 r(2 + A) 

dl ' cf (47r)^/2 r(2 + A) ' 

where the foUowing coefficients are used: 

2 1 - A 



C2 = 



2 



(27) 



(28) 



(477)^^/2 r(2 + A) ' 
~ _ 2 1 

^3 - - r(l + A) ' (29) 

d\ = -C3 , 

e 4 1 



(30) 



Ci (477)^^/2 r(l + A) ■ 

For the usual short range interactions A = 0, 

dl TT 

— = -/?(v) = vW , 

dl TT 

where 

ci = 1/8, C2 = I/tt, C3 = -I/tt, dl = I/tt (31) 

were used, e = ee is the order of e for the short range interaction and is dropped in the 
RG equation. 

The RG equations (28) for A tell us that the disorder is relevant for d> dc = 2+2A, 
irrelevant for d < dc and marginally irrelevant for d = dc- For dc < d < 3 + A, the pure 
fixed point W* = is unstable. However, there is a stable fixed point which is given 
by W* = e/Q where Q = (47^)^/2 r(2+A) • -^^^"^ we can read off the dynamical 
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exponent z = 1 + (4^)d/2 r(2+A) ^*' ^^^^ specifically, z = 1 + e/2 at this stable fixed 
point. For d < dc, the disorder is irrelevant and the pure fixed point is stable. Therefore, 
we can expect the direct superfluid-Mott insulator transition for d < dc- 

Now let us look at the short range interaction case A = 0. From Eq. (30), we can see 
that the disorder is relevant for cZ > 2 and irrelevant for d < 2. Therefore, there is only 
the unstable pure fixed point and the RG flow goes to the strong disorder regime for d = 2 
and slightly larger than two. For d <2, the pure fixed point becomes stable. There is also 
an unstable fixed point which is given by W* — |e|7r/128. The dynamical exponent at the 
unstable fixed point is ^ = 1 + |e|/2. 

III. RENORMALIZATION GROUP ANALYSIS FOR 
THE COMMENSURATE DIRTY BOSON MODEL 

Now we are going to study the generic model (3) with u = w = Here, we consider 
the usual short range interaction A = 0. The disorder characterized by 'u(x) in Eq.(3) 
will be considered in addition to the w{yi) type disorder. This means that wc have to 
consider more diagrams that are generated by this new disorder. The additional self 
energy correction due to ■u(x) type disorder is given by the diagram in Fig. 4 

d'^k 1 



J {27ry up- + fo(g - ky 



^6 

_ _r(i -rf/2) 

^4 - (47r)d/2 ' 
The new bare two point vertex F^^^^^ up to {j^Y' oi'der is 

r£L = + ^0?' - Si(9,a;) - S2(9,a;) 
.2 



= ^C4a;2a;'^-2 , (32) 



d-2 



(33) 



,22 (^_^^C22- 

-hfog |i 2 „2^ 



It is clear that we need to do e = d — 2 expansion in order to handle the divergences. Let us 
introduce additional dimensionless measures of the disorder Un = ^A'^"^ and U = ^Li'^~''^ 
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in the bare and the renormahzed theory. Adding appropriate counter terms and using (5) 
(6) (7), we can again obtain the following equations for Z and v 



C2 C3 



(34) 



Wv + c^Uv . 



There can be 5 additional diagrams that contribute to T^^^^. Let us identify T^^^^ ^, 
rw,6' ^hLj^ ^hL,8^ ^hL,9 as the diagrams in Fig.5 (a), (b), (c), (d), (e) respectively. 
First of all, let us evaluate the new vertex V2 of Fig. 6 (a) 



V2 = -Uouj' 
Uo 



d'^p 1 1 



d-A 



d-6 



UJ 



1 r(2-(i/2) r2(d/2- 1) 
3 - (4^)d/2 r{d-2) ' 

1 r(2 - d/2) T{d/2 - 2) V{d/2 - 1) 



d 

ci4 = 



(47r)^/2 V{d - 4) 

where fog S> a; is assumed. The bubble lis of Fig. 6 (b) is given by 



n3(g,a; = 0) = -[/o 



d'^k d'^p dv 



1 



{2ttY (27r)rf 2tt v^{k - q)^ + v^k^ + 
1 1 



X 



VQ{k — p)"^ + v"^ VQ{k — p — q^ + v"^ 



(35) 



(36) 



where 



^l-d/2(i__,)l-d/2 



2 (477)^^+1/2 'J^ [x{l-x)z + y{l-y){l-z)]'/^-d 



(37) 



is a convergent integral. We can see that the bubble Us is smaller by the factor e ( putting 
d = 2 + e ) than Hi and 112 so that the contribution from T^^^.^ g is higher order in e and 
we can neglect it. 
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Using the calculated V2, we can evaluate F^l^ as 



.(4) 
bare 



/ V bare,i 



1 



,i + ^o(ni + n2 + n3) 



Wo + 2{Vi + V2) 



1 



1 + ^oHi 



Wo 



goiUi + U2] 

p(4) , p(4) 
bare,w bare,u ' 



Wo + 2(Vi + V2) 



90^1 



VFo-?7oa;^(l + 2Fi + 2F2) 



where 



(4) _ ^^^2{3-d) 
^ bare,w _2„2 ^ ^ 



4g 



d-4 



+ (i4g' 



d-6 



vogi 4 Z2 

,4- 



+ 2 



Wovi di 2-d 



2 2 
9o q 



a; 



r(4) 

bare,M 



-Uoto' 



1-2 



d-4 



d-6 



+ 2 



WoVq di 2-d 

2 2 ^ 



where (^2 = for the short range interaction is used. 

We take the foUwoing renormalization condition for T^'^^ = T^^ + 



r(4) 

w 



wf_ 

22/^ 



2(3-d) 



gi=g4=M, 92=93=0, linia^o t^l='^2=u;3='i^4=a^'A' 



= - u . 



(38) 



(39) 



(40) 



Adding appropriate counter terms and using the renomalization condition, we obtain the 
foUwoing equations 



d-2 



1 + 2^ di] In - 



d-2 



W e ~ A ~ A 

1 + ^ (C3 + C2 + 2— - 2di) In - - t/ C4 In 



t/o UZ- 



-'1 



2-d 



2-d 



Cl 



1 - 2tl^^ In ^ 



(41) 



W 



1 - ^ C3 + C2 + 2c/i In - + [/ C4 In - 



A 



A 



1^ 
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Using ij, independence of the bare parameters, we obtain the following equations for W 
and U 

^-^W = P(W) = -eW +\ C2 + £3 + 2— - 2di - C4 UW , 

li^U = (3{U) = et/ - 4(c3 + C2 + 2dr)WU + C4 , 
where £4 = ec4 and . For small e, the coefficients of these equations are given by 

Cl = 1/8, C2 = I/tt, C3 = -I/tt, C4 = l/(27r), di = I/tt . (43) 

We can ignore e which is the order of e. Now the RG equations for W, U and v are given 

by 

^ = -P{W) = eW+ ^UW + ^-^W' , 
at 2n TT 

^ = = -6?7 + —WU - , (44) 

dl TT 27r 

/Qf \ 64 ~ 1 ~ 
— = -P{v) = vW - —vU , 

at TT ZTT 

where Z is agian the logarithmic measure of the RG flow. The RG flow is drawn in Fig.7 
(a) for e > 0, (b) for e < 0. 

For (i < 2 (e < 0), we can see that there are three fixed points which are given by 
{W*,U*) = (0,0), (0, 27r|e|), (^|e|,0). The dynamical exponents at these fixed points 
are given by2;=l, l + |e|, l + ^|e| respectively. Note that all of these fixed points are 
essentially unstable because they have at least one relevant direction in the W — U plane. 
Among these three fixed points, {W*, C/*) = (0, 27r|e|) is a non-trivial fixed point that has 
the least number (one) of relevant directions. However, we need to fine tune the strength 
of the w(x) type disorder in order to reach this fixed point. For 2 < d < 3 {e > 0), there is 
only the trivial fixed point {W*, U*) = (0, 0). At this fixed point, U is relevant and W is 
irrelevant so that the RG flow goes to the strong randomness regime where our RG scheme 
breaks down. 

Therefore, in both cases, there is no stable non-trivial random fixed point near the 
pure XY fixed point. This means that the SF-insulator transition should be discribed by 
possible strong random fixed points. The conclusion we can deduce from these results is 
that, near d = 2, the direct SF-MI transition is unlikely to happen in the whole w — u 
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plane even for the weak disorder. This result essentially support the original picture of 
Fisher et a/.[l,2] (Fig.l (a)) that the SF-insulator transition should alwalys occur from the 
Bose glass phase rather than from the Mott insulator even for the commensurate case. 



IV. 1/N CORRECTION 

In this section, we are going to investigate the effects of the 1/N correction to the 
RG equation. The 1 /N correction due to the 0^ interaction in the previous commensurate 
dirty boson model will be considered. The most important 1 /N correction that can affect 
the RG equation enters in the coefficient of the W term in Eq.(44). This 1/N correction 
can be read off from the scaling dimensions of (f)^(f) at the pure fixed point [20]: 



= 2 - 77 , (45) 
where 

As a result, the critical dimension for W is changed from 2 to 2 — 27/. If the convention of 
e = d — 2 was taken, the RG equation becomes 

dl 2tt tt 

§ = .S+'-^WV-^U\ (47) 

dl TT 2n 

dv 64 ~ 1 ~ 
— = vW vU . 

dl TT ZTT 

Note that the density p has no anormalous dimension [16] and the RG equation for U is 
not affected at this order. 

Now there can be three possible RG flows that depend on the dimensionality of the 
system and N. For e = d — 2 > 0, there is only one fixed point which is the unstable 
trivial fixed point {W*, U*) = (0, 0). The RG fiow is given by Fig.7 (a) and it flows to 
the strong disorder regime. For e < and |e| < 2r/, there can be two fixed points which 
are {W*,U*) — (0,0), (0, 27r|e|) with the dynamical exponents z — 1, 1 + |e|. But all of 
them are unstable and the RG fiow is given by Fig. 8. On the other hand, for e < and 
|e| > 2?7, there is one more fixed point which is given by {W*, U*) = ( ifg(|e| — 2?7), 0) with 
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2; = 1 + |(|e| — 2r]). All of the three fixed points are still unstable and the RG flow goes to 
the strong disorder regime. The RG flow is again given by Fig. 7 (b). Therefore, after the 
inclusion of the correction, there is no stable weak random fixed point. The direct 
SF-MI transition is unlikely in the whole w — u plane. 

V. SUMMARY AND CONCLUSION 

We study a large N generalization of the commensurate dirty boson problem. The 
expansion allows us to treat interaction effects properly. On the other hand, the disorder is 
assumed to be weak and the pertubation in the strength of the randomness is performed. 
In order to understand the behaviors of bosons in this model, we need two types of disorder, 
ie., the random coefficient of the density term cj)^ 5o0) ■J^(x), and the random coeffient of 
the quadratic term (f)^(f)^ w{yi). 

For a restricted model with i((x) =0 which has an additional particle-hole symmetry, 
we introduce more general long range interactions V{q) = go/q^, (0 < A < 1/2) ( which is 
not a genuine long range interaction as mentioned in Sec. II ). The critical dimension for 
w{x) type disorder is found to be = 2 + 2A and we performed e = d — 2 — 2X expansion. 
It is found that A = case and A 7^ case show different behaviors and they are not 
continuously related in the e = d — dc expansion. 

First, let us look at the case of A 0. For d > dc, the pure fixed point is unstable, but 
there is a stable fixed point which governs the transition. For d < dc, the pure fixed point 
is stable and the disorder is irrelevant. Therefore, we expect the direct SF-MI transition. 
For d < dc, the critical point is the same as that for the pure system. For d > dc, SF-MI 
transition is described by a new non-trivial fixed point. 

For the short range interaction (A = 0), the disorder is relevant for d>2 and irrelevant 
for d < 2. There is only the unstable pure fixed point and the RG flow goes to the strong 
disorder regime for 2 < d < 3. For d < 2, the pure fixed point becomes stable and there 
is also an unstable non-trivial fixed point. Therefore, we expect a direct SF-MI transition 
for d < 2 but none for 2 < (i < 3. 

We consider the usual short range interaction for the general case of it(x) but 
u = (the commensurate dirty boson problem). The critical dimension of both types of 
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disorder is found to be d = 2 (at (1/Ar)° order) and we performed e = d — 2 expansion 
{d < 3). For d < 2, we have three fixed points and they have at least one relevant direction 
in the W — U plane. In the case of 2 < d < 3, the RG flow is governed by the w{:k.) 
type disorder and it flows to the strong disorder regime which cannot be reached by the 
pertubation in the strength of the randomness. Therefore, we expect that the direct SF-MI 
transition is unlikely to happen near two dimensions. 

The effects of the 1/N correction are considered. For the commensurate dirty boson 
model, the correction due to the (p"^ interaction does not change the qualitative features 
of the problem. There is still no stable weak random fixed point for d ^ 2 and the direct 
SF-MI transition is unlikely. 
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FIGURE CAPTIONS 

Fig.l Three possible phase diagrams for the dirty boson system described by Eqs. (2) and 
(3). Our RG results favor the case (a). The boson density commensurates with 
the lattice in the MI phase and on the dotted line. In (a) there is no direct SF-MI 
transition. In (b) direct SF-MI transition is described by a tricritical point. Both of 
the BG and the MI phases are insulators, but the former is gapless and the latter has 
a finite gap. The commensurate dirty boson model is defined on the dotted line and 
its extension in the MI phase. 

Fig. 2 (a) Self energy correction due to the w{x.) type disorder. The wavy line represents the 
interaction and the dashed line with a cross x means the impurity average. W means 
w{'x.) type disorder, (b) The polarization bubble Hi. (c) The renormalized interaction 
( thick line ) due to the w(x) type disorder. 

Fig. 3 (a), (b), (c), (d) are diagrams that contribute to the four point function up to W"^ 
order, (e) The renormalzed polariztion bubble Hi. The double line means the renor- 
malized full propagator, (f) The polariztion bubble 112. (g) The vertex Vi. 

Fig. 4 The self energy correction due to the •u(x) type disorder. The dashed line with a cross 
X means the impurity average. U means the tt(x) type disorder. 

Fig. 5 Additional diagrams that contribute to the four point function in the generic com- 
mensurate dirty boson model up to 2nd order in the impurity stregth. 

Fig. 6 (a) The vertex V2. (b) The polarization bubble lis. 

Fig.7 RG flows (a) for 2 < d < 3 and (b) for d < 2. 

Fig. 8 RG flows which has the 1/N correction for d < 2 and |e| < 2r]. 
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